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Abstract: Global fits to charmless B —>■ PP decays in the framework of fiavor 
SU(3) symmetry are updated and improved without reference to the sin 2/5 mea- 
sured from the charmonium decay modes. Fit results directly constrain the (p, f/) 
vertex of the unitarity triangle, and are used to predict the branching ratios and CP 
asymmetries of all decay modes, including those of the Bs system. Different schemes 
of SU(3) breaking in decay amplitude sizes are analyzed. The major breaking effect 
between strangeness-conserving and strangeness-changing decays can be accounted 
for by including a ratio of decay constants in tree and color- suppressed amplitudes. 
The possibility of having a new physics contribution to Ktt decays is also examined 
from the data fitting point of view. 
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1. Introduction 

In the standard model (SM) of particle physics, CP violation in the quark sector is 
postulated to be purely derived from the Kobayashi-Maskawa mechanism, in which a 
3x3 unitary CKM matrix Vckm with three mixing angles and one CP- violating phase 
is used to describe charge-current weak transitions between the up-type and down- 
type quarks [Q]. The unitarity condition on the first and third columns of the CKM 
matrix is often used to form a triangle on the complex plane because the lengths of 
its three sides are of the same order. An important program in flavor physics is to 
constrain this unitarity triangle using as many independent experimental inputs as 
possible, for both determining standard model (SM) parameters with high precisions 
and discovering any possible new physics effect. A lot of progress has been done in 
this direction with the help of a huge amount of B meson data collected in the past 
few years at the S-factories 0, ^. 

Although charmless modes are rare processes in B decays, they are very sen- 
sitive to the smallest CKM matrix elements Vub and Vu through decay amplitudes 
and mixing, respectively. Moreover, they provide us information about the weak 
phases associated with these two matrix elements. Some theoretical analyses have 
been done in recent years to globally fit to charmless B PP and VP decay data 
in the framework of QCD factorization 0] and fiavor SU(3) symmetry 0, ^, ^ §]. 
Here P and V refer to pseudoscalar and vector mesons, respectively. Since the weak 
phase (3 (0i) is more accurately determined from the time-dependent CP asymmetry 
analysis of Bd — > {cc)Ks modes, the result is usually used as an input in the theory 
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parameters. With more modes being observed and measured at higher precisions, 
it becomes possible to use purely rare decays to provide a completely independent 
determination of the unitarity triangle vertex (p, fj), expressed in terms of the Wolfen- 
stein parameters 0, without reference to the charmonium modes. It is therefore one 
objective of the current analysis to see whether the charmless B decay data alone 
also provide a CKM picture consistent with other constraints. 

In this paper, we perform fits to the available charmless B —>■ PP decays us- 
ing the flavor diagram approach . The fitting parameters include the Wolfenstein 
parameters A, p, and f/, magnitudes of different flavor amplitudes, and their associ- 
ated strong phases. To take into account SU(3) breaking, we also include breaking 
factors of amplitude sizes as our fltting parameters in some flts. Generally speaking, 
our flts render an area of the (p, fj) vertex slightly deviated from but still consistent 
with that obtained from other constraints. Aside from a decay constant ratio, the 
SU(3) breaking is seen at O(10)% level. From the extracted ranges of theory param- 
eters, we predict the branching ratios and CP asymmetries of all decays using flavor 
SU(3) symmetry, including the Bs system. The latter will be compared with data 
already or to be measured at the Tevatron, large hadron collider (LHC), and KEKB 
upgraded for running at T(55'). 

This paper is organized as follows. In Section H, we introduce the notations used 
in this analysis, including the fltting parameters. Flavor amplitude decompositions 
of the rare decay modes, along with the available experimental data on branching 
ratios and CP asymmetries, are summarized in this section. The fltting schemes and 
results are presented in Section |^, where predictions and outlook of yet-seen modes 
are also given. Finally, Section ^ summarizes our flndings. 

2. Formalism and Notation 

The magnitude of invariant decay amplitude A for a decay process B — > Mi M2 is 
related to its partial width via the following relation: 

T{B^M,M,) = -^\AW (2.1) 

where p is the 3-momentum of the flnal state particles in the rest frame of the 
B meson. To relate partial widths to branching ratios, we use the world-average 
lifetimes r+ = (1.638 ± 0.011) ps, r° = (1.530 ± 0.009) ps and tb, = (1.466 ± 0.059) 
ps computed by the Heavy Flavor Averaging Group (HFAG) fl^. Unless otherwise 
indicated, for each branching ratio quoted we imply the average of a process and its 
CP-conjugate one. 

To perform the flavor amplitude decomposition, we use the following quark con- 
tent and phase conventions for mesons: 

• Bottom mesons: B^ = db, = bd, B^ = ub, B^ = —bu, Bg = sb, B^ = bs; 
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• Pseudoscalar mesons: 7r+ = ud, 7r° = [dd — uu)/\/2, vr = — du, = us, 
= ds, K = sd, = —su, rj = {ss — uu — dd)/y/3, rj' = {uu + dd+2ss) /\/Q; 

The f] and t]' mesons correspond to octet-singlet mixtures 

?7 = cos 6*0 - r/i sin 6*0 , (2.2) 
= ?78 sin 6^0 + r^i cos 6'o . (2.3) 

As shown in Ref. 0, varying the mixing angle 6q does not improve the quality of 
fits. Therefore, here we fix Oq = sin~^(l/3) ~ 19.5° according to the above-mentioned 
quark contents of t] and rj'. 

We list flavor amplitude decompositions and averaged experimental data for 
B — s> PP decays in Tables |I| and |^. Values of measured observables are obtained 



by weighted-averaging over the results of the BaBar |jT5|, 0, |T7|, |T8|, |l^, |20|, ^ 



|23|, Belle [||, |2|, |6|, g^, |2|, H, CLEO |0|, |T|, |3|, |3|, g, and CDF |3|, ^ 



37[1 Collaborations. The standard deviation is scaled by the scale factor S (whose 
definition can be found, for example, in Ref. [^) if it is greater than 1 in order 
to take into account the discrepancy among different experimental groups. These 
include: Br{TT-7r^) {S = 1.3), Br^-R-T]) {S = 1.1), Br{7i-r]') {S = 1.4), ^(tt+tt-) 
{S = 2.6), 5r(7rV) (S = 1.4), AcpiK-n^) {S = 1.3), Br{K-7]) {S = 1.3), Br{7]'K^) 
{S = 1.3), and A{ri'Ks) {S = 1.4). Amplitudes such as annihilation and exchange 
diagrams are expected to be small and therefore neglected in the calculation. 

In the present approximation, we consider five dominant types of independent 
amplitudes: a "tree" contribution T; a "color-suppressed" contribution C; a "QCD 
penguin" contribution P; a "flavor-singlet" contribution S, and an "electroweak 
(EW) penguin" contribution Pew- The former four types are considered as the 
leading-order amplitudes, while the last one is higher order in weak interactions. 
There are also other types of amplitudes, such as the "color-suppressed EW penguin" 
diagram P^w^ "exchange" diagram E, "annihilation" diagram A, and "penguin anni- 
hilation" diagram PA. Due to dynamical suppressions, these amplitudes are ignored 
in the analysis. This agrees with the recent observation of the B^ — ^ K^K~ decay. 

The QCD penguin amplitude in fact contains three components (apart from the 
CKM factors): Pt, Pc, and Pu, with the subscript denoting which quark is running in 
the loop. After imposing the unitarity condition, we are left with two components: 
Ptc = Pt — Pc and Ptu = Pt ~ Pu, integrating out the t quark from the theory. 
For simplicity, we will assume the t-penguin dominance, so that Ptc = Ptu and are 
denoted by a single symbol P. The same comment applies to other penguin-type 
amplitudes {e.g., Pew and Pew) ^d^- 

In physical processes, the above-mentioned flavor amplitudes always appear in 
specific combinations. To simplify the notations, we therefore define the following 
unprimed and primed symbols for AS" = and \AS\ = 1 transitions, respectively: 

J. — vurri (vu I -yc \ pC J./ — V'"^ T fV"- _L \ pC' 

^ — -'dfe-^ \^ db ~^ ^ dbJ-'^EW ^ — -'sfeSt-i yi I gfjjjr^^r , 
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Table 1: Flavor amplitude decompositions for strangeness-conserving B PP decays. 
Measured branching ratios (in units of 10~^) and CP asymmetries are given in the last two 
columns. For those modes with time-dependent CP asymmetries, A and S are listed in 
the first and second rows, respectively. 

p^-(yi + Fj,)fp-ipfJ , p'^-{y;, + y,: 




s^-(F,UFl)(^S--Pwj , s'^~{Y,l + Y^ 

e ^ Y-E - {Yl + Y^,)PA , e' = Y^E - (F,", + F/JPA , 

where Fj, = VqiqV*,^. Unless they are leading contributions, amplitudes such as e 
and pa are omitted from Tables and 0. 

One differnce between the current analysis and our previous work |^ is that the 
CKM matrix elements associated with the amplitudes are factored out here. The 
strong phases, however, are still absorbed in the amplitudes. Notice that when going 
from AS* = to \^S\ = 1 transitions, we explicitly include SU(3) breaking factors 
^t,^c, and for the T, C, and P amplitudes, respectively. In the naive factorization 
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Table 2: Flavor amplitude decompositions for strangeness-changing B — > PP decays. 
Measured branching ratios (in units of 10~®) and CP asymmetries are given in the last two 
columns. For those modes with time-dependent CP asymmetries, A and 5 are listed in 
the first and second rows, respectively. 

approximation, these SU(3) breaking factors are all equal to ^ = fx/ f-K = 1-223 
[PH . As an example, using the above-defined notations we have for the 5° K^tt^ 
decay: 

This can be obtained from the complete set of flavor amplitude decompositions given 
in Table ^. 

The CKM factors used in the analysis are given in terms of the Wolfenstein 
parameterization of the CKM matrix to O(A^). Since A has been determined from 
kaon decays to a high accuracy, we will simply use the central value 0.2272 quoted 
by the CKMfitter group as a theory input, and leave A, p = p(l — A^/2), and 
fj = ri{l — A^/2) as fitting parameters to be determined by data. 

A relation between the sizes of EW penguin amplitude and tree-type amplitudes 
has been found in Ref. [0] where the Fierz transformation is used to relate EW 
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penguin operators with tree-level operators. Explicitly, 



Pew = -SEw\T + C\e' 



EW 



(2.5) 



where Spj^^ is the strong phase associated with Pew- In the SM, 




0.0135 ±0.0012 . 



(2.6) 



In our fit, we will leave it as a free parameter to test how well the above relation 
holds. 

For the B meson decaying into a CP eigenstate fcp, the time-dependent CP 
asymmetry is written as 



where Am^ is the mass difference between the two mass eigenstates of B mesons 
and t is the decay time measured from the tagged B meson. 

3. Fitting Analysis 

To see the effects of SU(3) symmetry breaking, we consider the following four fitting 
schemes in our analysis: 

1. exact flavor SU(3) symmetry for all amplitudes {i.e., = C,c = = 1); 

2. including the factor fx/U for |r| only {i.e., = fx/U and ^^ = ^p = 1); 

3. including the factor fx/f-K for both |T| and \C\ {i.e., 6 = = fx/fn and 

= 1); and 

4. including a universal SU(3) breaking factor ^ for all amplitudes on top of 
Scheme 3 {i.e., ^t = ^^ = ^fx/U and = 0- 

To reduce the number of parameters, we assume exact flavor SU(3) symmetry for 
the strong phases in these flts. As a phase convention we set the tree amplitude to 
be real and positive, i.e., 6t = 0. 

In addition to the observables in B —>■ PP modes, we also include \Vub\ = 
(4.26 ±0.36) X 10~3 and = (41.63 ±0.65) x lO^^ as our fltting observables. Here 
we take the averages of their values measured from inclusive and exclusive decays as 
quoted in Ref. 0. They mainly help flxing the values of A and a/pM-^. We will 
discuss below how our flt results and predictions change should we use a lower value 
of Vub in the following numerical analysis. 



Acp{t) 



r(B° ^ fcp) - r(i?° ^ fcp) 

T{B^^ fcp) + T{B<^^ fcp) 

S sm{AmB ■ t) + A cos{AmB ■ t) , 



(2.7) 
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3.1 Fits to modes with only vr, K mesons in the final state 

We start by fitting to the branching ratios and CP asymmetries of the tttt, tcK, and 
KK modes of B meson decays. This part of the analysis avoid uncertainties in the 
wave functions and singlet amplitudes associated with the rj and r]' mesons. 

Currently, there are 20 experimental observables. Along with \ Vub\ and \ Vcb\ men- 
tiones above, we have totally 22 data points. The number of theoretical parameters 
is 10 for Schemes 1 to 3 and 11 for Scheme 4. The best-fitted values of the parameters 
in their 1 a ranges along with the minimal values, Xmm' ^"^^ the different schemes 
are listed in Table 0. 



Parameter 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 


P 


n 1 Qq+0.042 

U.iOJ_Q Q3y 




n 1 Q4+0.041 


n 1 qq+0.039 
U.±00_o 035 


V 


0.401 ±0.030 


0.403 ±0.031 


0.404 ±0.031 


0.399 ±0.031 


A 


0.807 ±0.013 


0.807 ±0.013 


0.807 ±0.013 


0.807 ±0.013 



ITI n c;7Q+0.055 n c,7c,+0.055 n e;7^+0.055 n c,Q9+0.056 

I J I U.O^-3_o,o47 U.O^O_o.047 U.O^^_o.047 U.OOZ_o.049 

|C| 0.371 ±0.050 0.364 ±0.050 0.364 ±0.049 0.372 ± 0.051 

5c -57.6 ±10.3 -55.9 ±10.7 -55.8 ± 10.2 -56.3 ±10.1 

\P\ 0.121 ±0.002 0.122 ±0.002 0.122 ±0.002 0.117 ±0.008 

Sp -22.7 ±4.0 -18.8 ±3.2 -19.3 ± 3.2 -IS.GtH 

IP I n nil +0-006 r, nil +0006 r,r,ii+o-oo5 nnii+o-oo4 

\r'EW\ U.Uii_o.oo3 U.Uii_o.oo3 U.Uii_o.oo3 U.Uii_o.003 

6p,^ 2.2^1:° -W.Ollli -15.1 ±39.9 

e 1 (fixed) 1 (fixed) 1 (fixed) IMt^^f^ 

Sew 0.013 ±0.006 0.013 ± 0.006 0.013 ± 0.005 0.013 ± 0.004 

xLJdof 18.9/12 18.0/12 16.4/12 16.1/11 



Table 3: Fit results of the parameters for the vrvr, ttK, and KK modes in Schemes 1 
through 4 defined in the text along with the minimal value. The amplitudes are given 
in units of 10^ eV. 



Generally speaking, we obtain fairly stable results for the parameters, except 
for some small variations in the strong phases among the SU(3) breaking schemes 
considered here. The fit quality is best in Scheme 3, suggesting that it is better 
to include the SU(3) breaking factor fx/ fn for the T and C amplitudes when going 
from the strangeness-conserving modes to the strangeness-changing modes. Scheme 4 
introduces an additional SU(3) breaking factor ^, which is found to be about 1.04. 
The small difference in xLin between Scheme 3 and Scheme 4 turns out to reduce the 
fitting quality from 17% down to 14%. 

It has been suggested as a direct test of flavor SU(3) symmetry by compar- 
ing the extracted amplitude magnitudes of K^K and B^ K^K with 
B^ — s> K^TT^ because all of them have the same single penguin amplitude contribut- 
ing to the decays, except for the only difference in the CKM factors. This is verifled 
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experimentally according to the current data. Taking the averaged invariant ampli- 
tudes of B° K^K and ^ R-R as \p\ and comparing it with \p'\ obtained 
from — s> K^T^^ , one obtains \p/p'\ — 0.23 ± 0.02 consistent with the expected 
ratio |Kd/Ks|- This justifies our choice of not including the factor fx/ fw for SU(3) 
breaking in the penguin amplitudes. We also find that Xmin becomes worse when 
the factor is imposed on the QCD penguin amplitude. Therefore, the data indicate 
that to a good approximation factorization works better for the color-allowed and 
color-suppressed tree amplitudes {i.e., T and C). 

We observe a large |C|/|T| ratio of about 0.63 in all these fits. This is different 
from the expectation of the usual perturbative calculation within the SM. A possible 



explanation is given in Ref. |]T3| where next-to-leading order corrections to the inter- 
action vertex are found to enhance the color-suppressed amplitude C for Kit decays 
by a factor of 2 to 3 while keeping other amplitudes more or less unchanged. More- 
over, there exists a large relative strong phase of about (—56 ± 10)° between C and 
T. Therefore, it can play an important role in CP asymmetries. The values of these 
parameters are largely driven by the large branching ratio of vr^vr^ and Acpiir^K') 
being quite different from Acp{7^~^K~). 

A strong phase of about —20° is associated with the penguin amplitude. This is 
primarily demanded by the CP asymmetries of the tt^tt^ and K^tt~ modes, both of 
which involve the combination of t^'^ and p^'\ 

In all our fits, the parameter 6ew is seen to be very stable and close to the 
value in Eq. ( p.6|) , which shows that the EW penguin amplitude has a size roughly 
agreeing with the SM expectation. This is partly due to the fact that the latest 
data are moving towards the SM estimates and a larger best fitted 7 is favored, 
modifying the T — P and T — Pew interferences and enhancing Br{'ir'^K^). It does 
not necessarily mean that the possibility of new physics has been ruled out. As it 
will be shown below (see Section p.3|) , an electroweak penguin-like new contribution 
with a different CP phase can dramatically improve the quality of fits. Besides, we 
find that Pew has a strong phase of about —10° relative to T and about 45° to C. 

Since the 1 a and 95%CL ranges for the (p, fj) vertex have unnoticeable change 
in the four fitting schemes defined above, we only present as a representative our 
preferred set. Scheme 3, in Fig. |^. This also shows the stability of the {p,f]) values 
against SU(3) breaking. Scheme 3 gives the following results for the weak phases a, 
P, and 7: 

a =(83+^)°, or 69° < a < 96° (95% CL) ; 
/3=(26±2)°, or 2r</5<31° (95% CL) ; (3.1) 
7 =(72+^)°, or 62°<7<8r (95% CL) . 

As shown in Fig. |1|, the determined ranges of (p, f]) for all the schemes in our 
fitting are slightly higher than those given by the latest CKMfitter and UTfit results 
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Figure 1: Constraints on the (/>, ff) vertex using B — > tttt, Xvr, and KK data in Scheme 3 
defined in the text. Contours correspond to 1 o" and 95% CL, respectively. The crosses 
refer to the 1 a range given by the latest CKMfitter (open circle) and UTfit (filled square) 
results using other methods ||2|, |3| as a comparison. 



obtained using other observables [|, |^ : 
CKMfitter: p = 0.207; 



-0.036 
-0.034 5 



7] = 0.341 ± 0.023 ; 



UTfit: p = 0.197 ± 0.031 , f] = 0.351 ± 0.020 . 



(3.2) 



These values are indicated by crosses with an open circle and a filled square in the 
figure, respectively. The difference is to a large extent caused by the large value 
of \Vub\ used in our fits. Therefore, we obtain slightly larger phases /? and 7 but a 
smaller a. 

We now briefly comment on the effects of using a smaller value of \Vub\ in the fits. 
If we take \Vub\ = (3.50 ± 0.18) x 10~^ extracted from unitarity angle measurements 
only Xmin is improved by 1.1. P reduces to around 21° and 7 increases to about 
75°, with a almost unaffected. The magnitudes of \T\ and \C\ both become slightly 
larger, but their ratio stays the same. The other parameters do not change much 
either. The same features are also observed in fits with all PP modes to be discussed 
in the next section. However, it should be emphasized that our analysis purposely 
avoid inputs other than the charmless decay modes unless necessary (such as A and 
\Vcb\ mentioned above). We therefore do not use this smaller \Vub\ value in our main 
analysis, for it relies quite a lot on charmed B decays. 
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Our predictions for the branching ratios and CP asymmetries for all the ^ 
TTTT, Kn, and KK modes based upon the extracted parameters in Table ^ are given 
in Table |. 



Observable 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 




5.4 ± 1.1 


5.4 ± 1.0 


5.3 ±1.0 


5.3 ± 1.1 




1.6 ± 0.4 


1.6 ± 0.4 


1.6 ±0.4 


1.5 ± 0.4 


Br{7T-7T^) 


5.3 ± 1.2 


5.4 ± 1.2 


5.4 ±1.2 


5.4 ± 1.3 


Br{7r+K-) 


20.2 ± 1.0 


20.1 ± 1.1 


20.1 ± 1.1 


20.3 ±4.3 


Br{TT°K^) 


9.9 ± 1.0 


9.9 ± 1.0 


10.0 ±0.9 


10.1 ±2.3 


Br{n~K^) 


23.0 ± 1.1 


23.1 ± 1.1 


23.1 ± 1.1 


23.4 ±4.8 


Br{TT°K-) 


12.0 ± 1.2 


12.1 ± 1.2 


12.0 ± 1.1 


12.2 ±2.5 


BriK+K-) 














Br{K^K^) 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.2 


Br{K-K^) 


1.1 ± 0.1 


1.1 ± 0.1 


1.1 ±0.1 


1.0 ± 0.2 




0.32 ±0.07 


0.27 ±0.06 


0.28 ±0.06 


0.26 ±0.06 




0.47 ±0.15 


0.49 ±0.15 


0.49 ±0.14 


0.50 ±0.14 


Acp{n-n^) 


-0.01 ± 0.04 


-0.02 ± 0.03 


-0.01 ± 0.03 


-0.01 ±0.03 


Acp{n+K-) 


-0.08 ±0.02 


-0.09 ±0.02 


-0.09 ±0.02 


-0.09 ±0.02 




-0.07 ±0.03 


-0.08 ±0.02 


-0.09 ±0.03 


-0.10 ±0.03 


Acp{7i-K^) 














Acp{tt'K-) 


0.00 ±0.03 


0.00 ±0.03 


0.01 ±0.04 


0.02 ±0.04 


Acp[J^ ) 


U 


U 


U 


U 


A{K^K'') 














Acp{K-K^) 














5(7r+7r-) 


-0.580 ±0.130 


-0.585 ±0.130 


-0.584 ±0.130 


-0.565 ±0.141 


5(7r07r0) 


0.814 ±0.109 


0.812 ±0.108 


0.810 ±0.106 


0.786 ±0.113 




0.851 ±0.042 


0.850 ±0.041 


0.861 ±0.041 


0.858 ± 0.042 




-0.000 ±0.014 


-0.000 ±0.014 


-0.000 ±0.014 


-0.000 ±0.015 



Table 4: Predicted branching ratios (in units of 10~^) and CP asymmetries for B — > 
vrvr, A'vr, and KK modes in different schemes. The observables with vanishing entries are 
predicted to be identically zero in our analysis. Experimentally measured quantities, if 
any, are already given in the last two columns of Tables |l] and ^ for comparison. 



It is seen that the values in the table are quite stable and generally in agreement 
with the measured numbers or upper bounds within the errors. The largest comes 
from 5(7r°fCs), which is entailed to be even larger than S(cc)Ks- The CP asymmetry 
of K~7c^ is found to be close to zero, giving the second largest contribution to x^- As 
we will see in Section |3.3| , these discrepancies can be significantly reduced if a new 
amplitude is introduced. 
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The central values of ^(tt+tt^) are slightly smaller than the measured one, but 
deviate from zero at more than 4 a level. The predicted ^(tt^tt^) are noticeably 
different from zero. This is seen as a result of the absence of a dominant amplitude 
in the decay. We also find a sizeable S 0.8±0.1, to be verified experimentally. 

At this point, it may be helpful to compare the predictions of other approaches. 
The recent next-to-leading order (NLO) calculations in QCD factorization (QCDF) 
show some enhancements from strong penguin corrections at 0{a1) However, 
their predictions for CP-averaged decay rates of the t:K modes still tend to be lower 
than the experimental data using the default parameter set, and the observed large 
negative AcpiT^^ K~) is difficult to understand. The hard spectator-scattering cor- 
rections have also been calculated to this order, and are found to possibly have a 
significant impact on the tree-dominated B ^ ttti decays |^0|]. It remains to be seen 
if a complete next-to-next-to-leading order (NNLO) calculation can lead to a better 
agreement with the data. 

In the perturbative QCD (pQCD) approach, NLO calculations including vertex 
corrections, quark loops and magnetic penguins suppress Acp{7^^K~) while keep 
A.cp{'^~^K~) large enough with the correct sign |T^, both in good agreement with the 
data. Nevertheless, the corrections to tttt modes are ineffective so that the predicted 
Br{TT^7c^) remains small. 

A recent comprehensive analysis in the same set of observables has been per- 
formed in the heavy quark limit of QCD and in the soft-coUinear effective theory 
(SCET) This approach generally involves more hadronic parameters without 

the help of symmetries. At the LO, they predict a larger Br{7i~^K~), an Acpiir^K') 
close to Acp{ti~^K~), and an A^ti^K'^) opposite in sign to the data. 

In all these approaches, the predicted S{7i'^Ks) are close to the one from S{J/ ipKs) 
or even larger [l^, This leaves room for new physics interpretations if it is 
further confirmed by data at a higher precision. 

Our predictions for the branching ratios and CP asymmetries for all the Bg 
TTTT, Kn, and KK modes based upon the extracted parameters in Table ^ are given in 
Table They serve as another good testing ground for the fiavor SU(3) symmetry. 
Among all observables of the Bg decays only the branching ratio of the Bg K^K^ 
mode is observed at CDF [PS|. This mode has the same flavor amplitude decompo- 
sition as the B^ tt^K~ mode. Therefore, our predictions in this case are close 
to those for 5° tt^K^, apart from a small difference due to such factors as the 
masses and decay widths. They are all smaller than the measured value. Since this 
is only observed for the first time with somewhat large errors, a more precise deter- 
mination will be very helpful. Besides, this mode is predicted to have nonzero CP 
asymmetries according to the fits. 

As mentioned above, a good fiavor SU(3) symmetry relation has been observed 
between B^ — > K'^K and B~ tt^K^. It is therefore natural to use the Bg — > 
K'^K decay as another test because it also involves a single p' amplitude. We 
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Observable 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 
















Sr(7r07rO) 














Br{n+K-) 


5.0 ±1.0 


5.0 ± 1.0 


5.0 ± 1.0 


5.0 ±1.0 


Br{n^K^) 


1.5 ±0.3 


1.5 ±0.3 


1.5 ±0.3 


1.4 ±0.3 


Br{K+K-) 


18.9 ± 1.0 


18.8 ± 1.0 


18.8 ± 1.0 


19.0 ±4.0 


Br{K^K°) 


20.0 ± 1.0 


20.2 ± 1.0 


20.1 ± 1.0 


20.4 ±4.2 






























Acpin+K-) 


0.32 ±0.07 


0.27 ±0.06 


0.28 ±0.06 


0.26 ±0.06 


Ain'^Ks) 


0.47 ±0.15 


0.49 ±0.15 


0.49 ±0.14 


0.50 ±0.14 


A{K+K-) 


-0.08 ±0.02 


-0.09 ±0.02 


-0.09 ±0.02 


-0.09 ±0.02 


A{K^K^) 














S{n+n~) 














5(7r07r0) 
















0.340 ±0.202 


0.365 ±0.194 


0.359 ±0.193 


0.308 ±0.201 




0.147 ±0.022 


0.199 ±0.028 


0.198 ±0.028 


0.211 ±0.035 




-0.043 ±0.004 


-0.044 ±0.004 


-0.044 ± 0.004 


-0.043 ± 0.004 



Table 5: Predicted branching ratios (in units of 10^^) and CP asymmetries for Bs decays 
in different schemes. The observables with vanishing entries are predicted to be identicahy 
zero in our analysis. 



predict its branching ratio to be around 2 x 10^^. Moreover, the time-dependent CP 
asymmetries A and S associated with this mode are predicted to be identically zero 
and about —0.044 ± 0.004, respectively. The Bg K^K and K^K~ decays have 
also been discussed in the literature to study their correlation with the B^ —>■ tt+tt^ 
decay and their time-dependent CP asymmetries for identifying new physics 



if they deviate from the SM predictions. 

The Bg —>■ K^-K^ and K^K~ modes have the same flavor amplitude decom- 
positions as the B^ -k^-k^ and tt^K" decays, respectively. Therefore, they are 
predicted to have sizeable CP asymmetries due to the interference between tree and 
penguin amplitudes. 

The same color-suppressed and penguin amplitudes contribute to both 5° 
71^71^ and Bs tt^K^ modes. Neither of them is dominant in the decay processes. 
Therefore, we expect large CP asymmetries in the latter mode as well. Moreover, de- 
termining the branching ratio of the latter may provide some insight for the observed 
large branching ratio of the former. 
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3.2 Fits to all B ^ PP modes 

We further carry out the analysis with the inclusion of modes with rj and rj' in the 
final state. In this case, there are totally 34 experimental observables. The best 
fitted parameters are listed in Table 0. 



Parameter 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 


P 




0.087^H^^ 


0.087^H^^ 


096+° °^^ 
u.uyu_Q 


V 


0.377 ±0.027 


0.378 ±0.028 


0.379 ±0.027 


0.370 ±0.027 


A 


0.809 ±0.012 


0.809 ±0.012 


0.809 ±0.012 


0.809 ±0.012 



|r| OMlt',f,t 0M2t',f,l OMOt',tl OM9t'otl 

\C\ 0.426 ±0.048 0.418 ± 0.048 0.415 ± 0.047 0.436 ± 0.049 

6c -72.5 ±7.3 -70.4 ±7.5 -70.0 ± 7.3 -68.3 ± 7.2 

\P\ 0.121 ±0.002 0.121 ±0.002 0.121 ±0.002 0.110 ±0.008 

5p -17.8 ±3.2 -16.0 ±2.8 -16.4 ± 2.8 -15.9 ±2.6 

\P,^\ O.OUt'oZl OMlt'oZl O.OUt'ofol OMSt'ofol 

5p..,. -h^.^tf.l -47.7+^H -58.1+?^-? -57.6+^^-^ 



\S\ 0.0481°:°°^ 0.047+°;°°^ 0.047+°;°°^ 0.042 ±0.004 

5s -48.3 ± 10.6 -44.8 ± 10.2 -44.2 ± 9.8 -42.9 ± 9.3 

i 1 (fixed) 1 (fixed) 1 (fixed) l-lOiH? 



5ew 


0.014 ±0.006 


0.013 ±0.005 


0.014 ±0.006 


0.015 ±0.006 


xVdof 


37.4/22 


34.8/22 


32.9/22 


30.6/21 



Table 6: Fit results of the parameters for all PP modes in Schemes 1 through 4 defined 
in the text along with the associated minimal values. The amplitudes are given in units 
of 10^ eV. 

To fit all the PP modes, we have to include at least the flavor singlet amplitude 
S, whose importance for explaining the large branching ratios of the rj'K modes has 
been noticed and discussed in Refs. 46, 47 1. This introduces two more theoretical 



parameters 15*1 and 6s, the strong phase associated with 5", than the fits in Section 
However, the fitting quality in these schemes is seen to be much worse than before. 
Among the modes with t]^'^ in the final state, Br(r]K~), Sjj'Ks^ Acp{7T~r]') have 
the largest contributions to x^- 

Comparing the fitting results in Table ^ with those in Table ^ we see that the 
strong phases suffer from larger fluctuations among all theoretical parameters. The 
values of both p and fj are decreased, but their precisions improved. This leaves a 
smaller region for the (p, f/) vertex, with a jS consistent with other observations and 
a somewhat larger 7. The parameters |T| and \C\ become slightly larger; but the 
ratio |C|/|T| ~ 0.65 remains about the same. The magnitudes of P and Pew (or 
6ew) are seen to be relatively stable in both limited and global fits. The parameter 
^ increases from 1.04 to 1.10. 
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The singlet amplitude has a magnitude about 3 to 4 times \Pew\ in our fits. 
Moreover, its strong phase is close to ^p^,^ and about —30° from P. It is this feature 
that produces interesting interference patterns among the different modes involving 
rj and f]' . 



0.8- 



0.6- 



0.2- 



-0. 



- 
- 














n 




1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.2 



0.4 



0.6 



0.8 



1.2 
P 



Figure 2: Constraints on the (p, f]) vertex using all the PP mode data in Scheme 3 defined 
in the text. Contours correspond to 1 o" and 95% CL, respectively. The crosses refer to 
the 1 a range given by the latest CKMfitter (open circle) and UTfit (filled square) results 
using other methods ^ as a comparison. 



Scheme 3 in this case gives the following results for the weak phases a, /3, and 



a = (80 ± 6)° , or 69° < a < 92° (95% CL) ; 
/? = (23 ± 2)° , or 20° < /3 < 27° (95% CL) ; 
7 = (77 ± 4)° , or 69° < 7 < 84° (95% CL) . 



(3.3) 



The predictions of the Bu^d decay observables are given in Table |^. Vanishing 
observables in our approach are omitted to avoid an oversized table. The following 
observables deviate the most from the current data: S{t[^Ks)i Brln^TT^), Br^K^r])., 
Briji^Tx^), Sijj'Ks), and Acpiji^rj) (hsted in the order of their contributions to x^)- 

As in the limited fits in Section |3.1| , the global fits also prefer sizeable CP asym- 
metries for the vr'^TT^ and vr+vr" modes. The branching ratios of the yet measured 7777, 
rjrj' , and rj'rj' are all consistent with the current upper bounds. Their corresponding 
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direct CP asymmetries are predicted to be large. However, measuring them will 
require more work. 

The current branching ratios of TT^rj and tt^t]' has a factor of 2 difference in the 
central values, though the errors are still large. Our results, on the other hand, show 
that they are equal to each other in all schemes. A^rjKs) has not been observed, but 
is 2 cr away from zero in our analysis. 

The predictions for the Bg modes are given in Table |[ Many of the discussions 



regarding the modes with tt and K mesons in the final state in Section |3.1| can be 
applied here. Thus, we only concentrate on the modes with t] and/or rj' in the final 
state. 

As in the cases of tj'Ks and B^ K't]', the constructive interference 

between p' and s' makes the B^ — ^ rj't]' the one with the largest branching ratio, 
~ 50 X 10~^, among all. The same effect is seen in the B^ 7777' decay as well. In 
contrast, a destructive effect occurs in the B^ — * rji] decay, so that its branching ratio 
is only ~ 2 x 10"^ 

The B^ i]'^s decay is another place where the constructive interference be- 
tween the QCD penguin and singlet penguin amplitudes plays an important role. 
Although small in magnitude for AS* = transitions, they can interfere with the 
color-suppressed amplitude to give potentially observable time-dependent CP asym- 
metries, both predicted at ~ 5a level. 

3.3 Fits with a new physics amplitude 

In expectation of possible new physics contributions to the Ktt decays to account for 
the observed branching ratio and CP violation pattern ^ |5l|, we try in 
Scheme 3 fits with a new amplitude added to these decays. More explicitly, a new 
amplitude = |A^| exp [i{(j)N + ^n)] is included in the B — > tt^K' and tt'^K^ decays 
in such a way that effectively, 

c' ^ y:,^,c - (y;, + r/j Pew + n . (3.4) 

This introduces three more parameters (|A^|, (pN, and 6n) into the fits. Here we 
assume that Pew is fixed relative to T + C through the SM relation. Xmin is found to 
decrease from 16.4 to 4.3 in the limited fit with only tt, K mesons in the final state. 
The new physics parameters are found to be 

|Ar| = i8+3eV, 07v = (92±4)°, and = (-14 ± 5)° . (3.5) 

The best fitted CKM parameters p, fj and A remain almost unchanged. The best 
fitted tree-type amplitudes have |T| = (0.55+°;°^) x 10^ eV, \C\ = (0.32+°;°^) x 10^ eV 
and 6c = (— 39;1^}3)°. The penguin amplitude \P\ is unchanged. 

Since both C and Pew have the same flavor topology, they always appear in pairs 
in c or c' for any physical decay process. Therefore, it seems difficult to determine 
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whether the new amphtude is associated with one or the other . Our results have 
\N\/\VcbV^s\ ^ 0.04 X WeY and |A^|/|KfeKs| ^ 2.2 x lOW, showing that the new 
amphtude is unexpectedly large. It is about five times bigger than |PEvy| or \C\. 
Since we assume that it only enters c' in the Kit modes instead of c in the vrvr modes, 
this result suggests that it behaves more like the electroweak penguin amplitude than 
the color-suppressed amplitude, for the former plays a much less important role than 
the latter in the strangeness-conserving modes. Although the latest data indicate 
only a mild deviation of t^K branching ratios from SM estimates, the large difference 
between Acp{i^^K~) and Acp{i^^K~) is unexpected, which require a much larger 
|C"|/|T'| than |C|/|T| in tttt modes [Q. Within the framework of SU(3) symmetry, 
a large new physics contribution is still possible. 

The above conclusion may look contradictory to what we have found in Sec- 
tion 3A, where |Pevk| is preferred by data to fall within the SM expectation, meaning 
that varying its value would not improve the fitting quality. But this is only because 
in the previous fit, the weak phase of Pew is fixed to the SM value. In the analysis 
of this section, the electroweak penguin-like new amplitude is allowed to have its 
own weak and strong phases. 

It should be emphasized that N is not added to modes with the contributing 
amplitude d in the global fits other than the Kit decays. It does not improve the 
minimal much to do so, for there is a pull between the Ki: and Kr]'^''' data such 
that |A^| is about a factor of 5 smaller than that quoted above. Therefore, it remains 
to be understood why the new amplitude does not help when modes with rj and r]' 
are taken into consideration as it should if fiavor SU(3) is respected. The solution to 
this question may rely on more precisely determined branching ratios of 7]K^'~ . 



4. Summary 

In this paper, we perform fits to the branching ratios and CP asymmetries of 
both limited and entire sets of the rare B PP decays. We consider the primary 
contributing fiavor amplitudes T, C, P, and Pew, each of which is associated with 
a distinct strong phase. The analysis is based upon fiavor SU(3) symmetry. We 
also include the fx/ fw ratio and an additional SU(3) breaking factor ^ to test the 
stability of our fits. 

One major result is the extraction of the vertex (p, f/) and thus the weak phases 
of the CKM unitarity triangle. This is complementary to other methods. The values 
of P and 7 obtained from our fits for modes without i]^'^ in the final state are generally 
larger than but still consistent within errors with those given by the overall fits of 
the CKMfitter and UTfit groups to other observables. Our fits to all the PP modes, 
however, result in a (3 similar to that given by both the CKMfitter and the UTfit 
groups but a larger 7. 
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The current PP data favor a large C with a strong phase of about —60° with 
respect to T. This is seen to be required by the large branching ratio of vr'^vr^ and 
the fact that Acp{Ti^K~) is different from Acp{'n:^K~). On the other hand, the 
size of electroweak penguin amplitude Pew is found to be consistent with the SM 
expectation. 

We also comment on the possibility of a new physics contribution to K-k decays. 
Our fitting analysis in this case prefers an electroweak penguin-like amplitude with 
sizeable magnitude and nontrivial weak phase given in Eq. ( ^.51 ). However, this 
amplitude does not respect flavor SU(3) symmetry. 

Using the parameters extracted from fitting, we make predictions for all the 
rare B —>■ PP processes. Moreover, we extend our predictions to the observables in 
Bs decays based on flavor SU(3) symmetry, whose experimental data will become 
available for comparison in the next couple of years from Tevatron Run II, LHCb 
and upgraded Belle experiments. 
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Observable 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 


-Br(7r+7r^) 


5.3 ±1.0 


5.3 ± 1.0 


5.3 ± 1.0 


5.3 ± 1.0 


Sr(7rV) 


1.7 ±0.3 


1.7 ±0.3 


1.7 ±0.3 


1.6 ±0.3 


Sr(7r-7r°) 


4.8 ± 1.0 


4.9 ± 1.0 


4.9 ± 1.0 


5.1 ± 1.1 




20.3 ± 1.0 


20.2 ± 1.0 


20.2 ± 1.0 


20.4 ± 4.3 


Brin^'K^) 


9.6 ± 1.0 


9.6 ±0.9 


9.6 ± 1.0 


9.8 ±2.3 


Br{TT-K^) 


22.6 ± 1.1 


22.7± 1.1 


22.7 ± 1.1 


23.1 ±4.8 


Br(Ti^K-) 


12.3 ± 1.2 


12.2± 1.1 


12.3 ± 1.2 


12.5 ± 2.7 


Br{K^K^) 


1.1 ±0.1 


1.1 ±0.1 


1.1 ±0.1 


0.9 ±0.1 


Br{K-K^) 


1.2 ±0.1 


1.2 ±0.1 


1.2 ±0.1 


1.0 ±0.2 


Br{'K^rj) 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.1 


0.8 ±0.1 


Br{TT^h/) 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.1 


0.8 ±0.1 


Br{7r~r]) 


4.6 ±0.6 


4.6 ±0.6 


4.6 ±0.6 


4.6 ±0.7 


Bri-K'T]') 


3.2 ±0.3 


3.2 ±0.3 


3.2 ±0.3 


3.0 ± 0.4 


Br{K^7]) 


1.4 ±0.2 


1.3 ±0.2 


1.4 ±0.2 


1.4 ±0.3 


Br{K^r]') 


65.3 ±5.2 


65.7 ±5.0 


65.5 ±4.8 


66.4 ± 13.0 


Br{K~Ti) 


1.5 ±0.3 


1.5 ±0.2 


1.5 ±0.3 


1.5 ±0.4 


Br{K-r]') 


69.2 ±5.5 


69.5 ±5.3 


69.3 ±5.1 


70.1 ±13.8 


BrijiT]) 


0.8 ±0.1 


0.8 ±0.1 


0.8 ±0.1 


0.8 ±0.1 


Br{ri'ri') 


0.4 ± 0.0 


0.4 ±0.0 


0.4 ±0.0 


0.4 ± 0.0 


Br{rjrj') 


1.2 ±0.1 


1.2 ±0.1 


1.2 ±0.1 


1.1 ±0.1 


^(Tr+Tr") 


0.27 ±0.06 


0.24 ±0.05 


0.25 ± 0.05 


0.22 ± 0.04 


^(7r°7r°) 


0.71 ±0.10 


0.70 ±0.10 


0.70 ±0.10 


0.67 ±0.09 




0.03 ±0.03 


0.02 ±0.03 


0.03 ±0.03 


0.04 ±0.03 


Acp{n+K-) 


-0.07 ±0.02 


-0.08 ±0.02 


-0.08 ±0.02 


-0.08 ±0.02 




-0.13 ±0.02 


-0.12 ±0.02 


-0.15 ±0.03 


-0.17 ±0.03 


Acp{7i~rj) 


-0.09 ±0.10 


-0.11 ±0.09 


-0.10 ±0.09 


-0.10 ±0.09 


Acp{n-r]') 


0.06 ±0.12 


0.04 ±0.12 


0.04 ±0.12 


0.02 ±0.11 


AivKs) 


0.13 ±0.07 


0.14 ±0.07 


0.16 ±0.08 


0.18 ±0.09 


Aiv'Ks) 


0.02 ± 0.00 


0.02 ± 0.00 


0.03 ±0.01 


0.03 ±0.01 


Acp{K-r)) 


-0.25 ±0.12 


-0.29 ±0.13 


-0.27 ±0.14 


-0.29 ±0.15 


A{r]r]) 


-0.77 ±0.11 


-0.78 ±0.11 


-0.76 ±0.11 


-0.73 ±0.11 


AivW) 


-0.55 ±0.13 


-0.55 ±0.13 


-0.55 ±0.13 


-0.58 ±0.13 


A{7]7]') 


-0.65 ±0.13 


-0.66 ±0.13 


-0.66 ±0.12 


-0.66 ±0.12 


S{7l+7l-) 


-0.533 ±0.135 


-0.533 ±0.135 


-0.532 ±0.134 


-0.513 ±0.138 




0.634 ±0.116 


0.655 ±0.111 


0.649 ±0.111 


0.614 ±0.118 




0.780 ± 0.041 


0.781 ± 0.041 


0.789 ± 0.041 


0.791 ± 0.041 


S{K^K^) 


-0.001 ± 0.031 


-0.001 ± 0.031 


-0.001 ± 0.030 


-0.000 ±0.017 


S{r)Ks) 


0.5 ± 0.06 


0.5 ±0.05 


0.45 ± 0.06 


0.39 ±0.07 


Sin'Ks) 


0.72 ±0.01 


0.72 ±0.01 


0.72 ±0.01 


0.71 ± 0.04 



Table 7: Predicted branehing ratios (in units of 10^^) and CP asymmetries for all the 
PP modes of B^'^ decays in different schemes. Vanishing observables in our approach are 
omitted. 
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Observable 


Scheme 1 


Scheme 2 


Scheme 3 


Scheme 4 


Sr(7r+7r~) 














Sr(7r07rO) 














Bri-K+K-) 


5.0 ±0.9 


5.0 ±0.9 


5.0 ±0.9 


5.0 ±0.9 


Bri-K^K^) 


1.6 ±0.3 


1.6 ±0.3 


1.6 ±0.3 


1.5 ±0.3 


Br{K+K-) 


18.9 ± 1.0 


18.9 ± 1.0 


18.9 ± 1.0 


19.1 ±4.0 


Br{K^K°) 


19.7 ± 1.0 


19.8 ± 1.0 


19.8 ± 1.0 


20.2 ±4.2 


Brln^f]) 








0.1 ±0.0 


0.1 ±0.0 


Sr(7rV) 


0.1 ±0.0 


0.1 ±0.0 


0.1 ±0.0 


0.1 ±0.1 


Br{K^r]) 


0.7 ±0.2 


0.7 ±0.2 


0.7 ±0.2 


0.7 ±0.2 


Br{Ky) 


3.3 ±0.3 


3.4 ±0.3 


3.4 ±0.3 


2.8 ±0.3 


Br{riri) 


2.0 ±0.4 


2.0 ±0.4 


2.0 ±0.4 


2.0 ±0.6 


Br{r]'rj') 


48.3 ±4.4 


48.6 ±4.3 


48.3 ±4.1 


48.9 ±9.8 


Br{r)r)') 


22.4 ± 1.5 


22.6 ± 1.4 


22.5 ± 1.4 


22.9 ±4.7 






























Acp{n+K-) 


0.27 ±0.06 


0.24 ±0.05 


0.25 ± 0.05 


0.22 ± 0.04 


A{n'Ks) 


0.71 ±0.10 


0.70 ±0.10 


0.70 ±0.10 


0.67 ±0.09 


A{K+K-) 


-0.07 ±0.02 


-0.08 ±0.02 


-0.08 ±0.02 


-0.08 ±0.02 


A{K^K^) 
















0.20 ±0.47 


0.32 ±0.48 


0.19 ±0.46 


0.18 ±0.45 




0.20 ±0.47 


0.32 ±0.48 


0.19 ±0.46 


0.18 ±0.45 


A{r^Ks) 


-0.24 ±0.13 


-0.27 ±0.13 


-0.26 ±0.12 


-0.22 ±0.11 


Aiv'Ks) 


-0.42 ±0.08 


-0.41 ±0.08 


-0.41 ±0.08 


-0.45 ±0.08 


A{vv) 


-0.20 ±0.03 


-0.20 ±0.03 


-0.24 ±0.04 


-0.27 ±0.05 


AivW) 


0.03 ±0.01 


0.03 ±0.01 


0.03 ±0.01 


0.04 ±0.01 


Aim') 


-0.02 ± 0.00 


-0.02 ±0.00 


-0.03 ±0.01 


-0.03 ±0.01 


«S(7r+7r-) 














5(7r°7r0) 














0/ 7"^ \ 


0.282 ± 0.158 


0.318 ± 0.153 


0.311 ± 0.153 


0.185 ± 0.167 


S{K+K') 


0.167 ±0.024 


0.217 ±0.030 


0.216 ±0.030 


0.244 ±0.037 




-0.041 ±0.004 


-0.041 ±0.004 


-0.041 ±0.004 


-0.040 ± 0.003 


SivKs) 


0.26 ±0.17 


0.24 ±0.17 


0.26 ±0.17 


0.26 ±0.16 


Sirj'Ks) 


0.40 ± 0.08 


0.39 ±0.08 


0.39 ±0.08 


0.45 ± 0.08 



Table 8: Predicted branching ratios (in units of 10 ^) and CP asymmetries for all the PP 
modes of Bg decays in different schemes. Observables with vanishing entries are omitted. 
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